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Abstract 

We study Penner type matrix models in relation with the Nekrasov partition 
function of four dimensional J\f = 2, SU{2) supersymmetric gauge theories with 
Np = 2, 3 and 4. By evaluating the resolvent using the loop equation for gen- 
eral (3, we explicitly construct the first half-genus correction to the free energy 
and demonstrate the result coincides with the corresponding Nekrasov partition 
function with general 0-background, including higher instanton contributions 
after modifying the relation of the Coulomb branch parameter with the filling 
fraction. Our approach complements the proof using the Selberg integrals di- 
rectly which is useful to find the contribution in the series of instanton numbers 
for a given deformation parameter. 
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1 Introduction 



In 2009, a remarkable relation, so called "AGT conjecture" |JJ was proposed between 
the Nekrasov partition function of A/" = 2 SU{2) gauge theory in 4D [2l [3] and the 
conformal block of the Liouville theory in 2D, relating the f2-background parameters 
ei and €2 and instanton expansion parameters with background charge of the Liouville 
theory and modular parameter of the conformal block. In addition, the vacuum expec- 
tation value of the Coulomb branch parameter (VEV of the adjoint scalar) is related 
with the momentum of the primary field of the intermediate channel. Soon after this 
relation was generalized to SU{N) gauge group in [1] and [HI E]. 

After this conjecture, various works have been performed using the Selberg integral 
and Jack polynomials such as in [H EJ |9l [TOl [TTl [12] and references therein. The 
Nekrasov partition function is obtained in the limit ei + €2 = corresponding to the 
c = 1 Liouville theory (which is also called /3 = 1 limit, which will be elaborated 
later). Nevertheless, the conjecture seems to go beyond this limit /3 7^ 1 and the 
Selberg integral provides a nice tool to this approach. 

Similarly related but a little different approach we are going to investigate in this 
paper is to view the partition function in terms of Penner type matrix model. This 
was initially proposed in [13] for four flavor case and generalized to less flavor cases 
in [H] noting that the Liouville conformal block can be reproduced in terms of the 
/3-deformation of hermitian matrix model. In this matrix model, the gauge parame- 
ters and the f2-background parameters determines the matrix couplings, the deformed 
parameter /3 and the size of the matrix A^. These matrix models and related topics 
have widely been studied in [T5]-[25]. 

The free energy of the /3-deformed matrix model is generally expanded in powers 
of the coupling g corresponding to the genus expansion 

00 

VlogZ = ^F„(A)^?", (LI) 

n=0 

where A is a parameter involved in the matrix potential and interpreted as a dynamical 
scale of the corresponding gauge theoryl^l When /3 = 1, the genus expansion terms F„ 
are vanishing for odd n. The planar free energy -Fo(A) was shown to be equivalent to 
the Seiberg-Witten prepotential [HI [TJ [21], and the genus one correction -Fi(A) was 
also shown to be consistent with the Nekrasov partition function llSj. 

When /3 7^ 1, the genus expansion terms Fn{A) in (11. ip is not vanishing also for 
odd n. This is easily seen in the loop equation of the /3-deformed one such as given 

■fWhen the corresponding gauge theory has four flavors, A should be identified with the exponential 
of a UV gauge coupling. 
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in [26]. In this paper, we apply the usual loop equation technique to the relatively 
simple model, namely, the /3-deformed matrix models for Np = 2,3,4 for SU{2) gauge 
group case and evaluate the half-genus correction -Fi(A), closely following the method 
employed in [21] but generalizing to (3 1 case and evaluate -Fi(A). In fact, the loop 
equation was initially studied in [23] and only an integral expression is presented. In 
this paper, we will calculate the explicit expression of the first half-genus correction 
to the free energy and compare the result with the corresponding Nekrasov partition 
function with general i7-background, including higher instanton contributions. 

The rest of this paper is organized as follows. In section |2l we briefly review 
the matrix models for Np = 2,3 and 4, and evaluate the resolvent of the matrix 
models using the loop equation for /3-deformed case. At the planar limit, the filling 
fraction is identified with the Coulomb branch parameter. However, it is noted that 
the relation is to be modified at the order of 0{g). In section [3l we concentrate on 
the Np = 2 case and evaluate the half-genus correction to the free energy. In section 
m and we generalize the previous argument to the Np = 3 and 4 cases, respectively. 
Section 6 is the summary and discussion. In appendix A, Penner type matrix model 
is constructed from the AGT conjecture. In appendix B, the derivation of the loop 
equation is reviewed for general (3. In appendix C, the explicit expression for the 
Nekrasov partition function is shown for Np = 2,3 and 4 for comparison with the 
matrix model results. 



2 Penner type models and half-genus corrections 

In this section, we describe the Penner type matrix models proposed in [131 E] to set 
up our approach. Then, we solve the loop equation and find the expressions of the 
resolvent for the planar and half-genus correction. In addition, the Coulomb branch 
parameter is given in terms of filling fraction. 



2.1 Penner type models for SU{2) gauge theories 



Penner type matrix model was proposed in [13] for the Np = 4 case and later gener- 
alized in [H] to Np = 2,3 cases. The partition function Z^atrix of the matrix model is 
defined as 



n/ ^^^j ^^^^^p 
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(2.1) 



where g is a coupling constant of the matrix model, /3 is a deformation parameter and 
'^N = n/<j('^/ ~ -^j) is the Vandermonde determinant. When (3 = 1, this reduces to 
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the usual hermitian matrix model with a coupling g and Xj is the eigenvalues of the 
hermitian matrix. 

According to the AGT conjecture, the parameters g and /3 are related to the Q- 
background parameters ei, €2 of the gauge theorjH 

ei = 2g^, £2 = (2.2) 

The matrix size is identified with the number of screening charges in Liouville theory 
and is related to the the mass parameters of the gauge theory. The details depend on 
the number of flavors as explained below. 

The free energy of the matrix model is defined by 

^matrix ^ 4^2 ^^^^^.^ ^ ^^^^^.^^ ^2.3) 

and is expanded in powers of g as 

^matrix ^matrix _j_ ^+ j-imatrix _j_ ^^g^^ ^2 4) 

where e+ = ei + €2 is the order of g and the half-genus correction F^'^^'^^^ is of our chief 
concern. In comparison to -Fi(A) in (11. ip . we see gFi = ^^matrix_ 

The explicit form of the potential V{z) depends on the number of flavors Np in 
gauge theory. When Np = 4, the potential is given by 

V{z) = ^mo + logz + mi log(z - 1) + m2 log(2 - g), (2.5) 

where q in (12.51) is identified with the exponential of a UV marginal coupling in the 
gauge theory. The mass parameters mo,mi and m2 (and with additional moo) are 
associated to the Cartan sub-algebra of 5*0(8) flavor symmetry and are related to the 
masses ni of the four anti-fundamental hypermultiplets by 

/ii = mi + moo, fi2 = mi-moo, /is = "^2 + mo, fi^ = 1712 - rrio. (2.6) 

We put potential term proportional to e+ in (12.51) so that the mass parameter 
relation in (12. 6p is maintained (see appendix |A]) . The matrix size A^ is determined by 
the relation 

^I^+^is + 2^/|3gN = 0. (2.7) 

which corresponds to the neutrality condition in the presence of the background charge 
of the Liouville theory. At the planar limit as explicitly shown in [T3] the e_|_ dependent 
correction can be neglected. 



^We use the notation g as one half of the topological string coupling gs 
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The case with Np = 3 is obtained if one takes the hmit ^ oo while keeping 
A3 = fi^q finite so that a single hypermultiplet is decoupled [271 E]- The resulting 
potential is given as (neglecting a divergent constant term) 

Viz) = (/i3 + y)log;2 + milog(^-l)-^. (2.8) 

Here A3 is a dimensionful parameter, and identified with the dynamical scale of the 
resulting gauge theory. The potential for Np = 2 is obtained from Np = 3 case further 
by taking the limit fi2 ^ 00 with (A2)^ = fi2-^3 fixed0 

V{z) = (/^3 + ^)logz + ^(^ + i). (2.9) 

Note that the relation (12.71) is not changed by the limiting process since fii and fis are 
required to be finite in this limit. 



2.2 loop equation and spectral curve 

We now discuss half-genus correction to the resolvent of the Penner type matrix mod- 
els. The resolvent is defined by 



W{z) ^ 7/3^(>^__), (2.10) 
which can be expanded in powers of g: 

00 

W{z) = Y,Wn{z)g^ (2.11) 



n=Q 



The resolvent satisfies the /3-deformed version of the loop equation [261 1211 129] . For 
our purpose of studying the half-genus correction, we may neglect 0{g'^) (see appendix 
E]) to get 

Wizf + ^W'iz) + Wiz)V\z)-^ = 0, (2.12) 
where f{z) is defined by 

It is noted that Wn{z) for odd n does not vanishing when /3 7^ 1, resulting in the 
"half- genus expansion", rather than the usual genus-expansion for the case of /3 = 1. 



^The sign of A2 is different from that in [M], but it is just the matter of convention. 
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To solve the loop equation fl2.12p . we divide the potential into two parts [2^: 

V{z) = V,{z) + ^Vl(z) (2.14) 
with Vq{z) and Vi{z) = log 2;, 0{g^). Putting the resolvent of the form 

Wiz) = Woiz) + '-±W^iz) + Oig'), (2.15) 



one has 



{WoY + V^Wo-^ = 0, (2.16) 



{2Wo + Vo)Wi + Wo + WoV; = 0, (2.17) 
whose solution is given as 



w.iz) ^ -yo(^)+^f^y+f(^\ (2.18) 

_ W^jz) + Woiz)V;iz) 



Now, the resolvent defines the spectral curve 

X = 2W{z) + V'{z) =xo + yxi + Oig^), (2.20) 

where xo is the the planar contribution 

(^0)^ = Kizr + f{z) (2.21) 

and Xi is the half-genus correction 

^ -x',iz) + V^'iz) + V^iz)Vliz) _ 

The above Penner type potentials present two cuts in the spectral curve and it is 
proposed in [13] that the Coulomb branch parameter a of the gauge theory is to be 
identified with the "filling fraction" of the matrix model 

a = — I xdz, (2.23) 

where A denotes a cycle surrounding a branch cut associated to x. This proposal 
is explicitly checked in [21j for /5 = 1 and fl2.23p reproduces a correct UV behavior 
a ~ ^/u where u denotes the VEV of the adjoint scalar in gauge theory. 

On the other hand, the half-genus corrections of the matrix model forces one to 
modify the identification into 

a = ^ (p xdz + ^. (2.24) 



2iTi 2 



This modification is necessary to keep the asymptotic behavior a ~ y/u in the UV hmit 
of the gauge theory. This is because the half-genus correction of the filhng fraction 
does not vanish as m — t- oo and the second term in fl2.24p cancels the non-vanishing 
half-genus correction. 

To elaborate on this modification, we redefine a = oq + + 0{g'^), so that 



ao = (p xodz, ai = (b xidz + 1. (2.25) 

2711 /_4 2tci 

Here ai can be put into more useful form if one uses (12.221) 

2711 27n /_4 xo{z) 

The first term gives just —1 due to the monodromy of the logarithmic function, which 
encircles a square- root branch cut associated with xq. This enforces one to modify 
the identification of the Coulomb branch parameter (I2.24p so that this contribution is 
canceled out. Thus, one has the expression 

a. = ^ / l?:(£)±i»£),, ,2.27) 

2711 /_4 xo{z) 

In the next sections, we evaluate the free energy of the matrix models by using oq 
and ai and compare with the Nekrasov partition function for Np = 2, 3, 4 cases. 



3 Nf = 2 case 

The potential (12. 9 p (omitting the subscript of A2 for simplicity) is given as 

Vo{z)=fislogz + ^(z + -], Vi{z) = logz. (3.1) 



2 V 

To find the spectral curve we evaluate f{z) (I2.13P in form [HI EI] 



f{z) = - + ^, (3.2) 
z z^ 



^ 'A 



where ci and C2 are given by 

ci = 2y^gNA = -{ft, + /i3)A, ^2 = 2^^5(7 ^ . (3.3) 

The exphcit form of ci is obtained due to the relation (12. 7p . Since ci has no half-genus 
correction, the planar spectral curve has the same expression as in the /3 = 1 case [21] : 

( ^2 , A 1 y ci C2 A^Piiz) 
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where P4,{z) is a polynomial of degree four of the form 

P^{z) = z + — \a^ + C2 - —\ z - —z + l. (3.5) 

This shows that the spectral curve is parameterized by a single complex paramter 
C2, which will be interpreted as the moduli parameter of the gauge theory [HI [2T] . 
Hereafter, we set /xi = /is = m just for simplicity. Then Pi{z) becomes 

P,{z) = z'-—z' + —[rn' + C2-—y-j^z + l. (3.6) 

The Coulomb branch parameter a is obtained from the filling fraction (I2.24p . The 
planar contribution was evaluated in [21] using the series expansion of hypergeo- 
metric function and has the form 

r- f m^ „ (^2 - em^A + 15m^) . , bi'iw? - Um^ A + 2lm^) . 

ao = V A 1 — - ^ -A-^ ^ ^A^ 

" V 4^2 64^4 256^6 

15(A^ - 28m2A3 + 294mM2 - 924mM + lOOlm^; g 



16384A8 



A« + c)(AiO) (3.7) 



where A = im? + C2 — The half-genus correction oi is easily obtained if one notices 
the relatioiJij] ai = — ^ | ^. ^^^^ because and ai have the form 



ao 



^£^V^c^. (3.8) 
a, = ^/ fi + .l^^^. (3.9) 



This follows from fl2.27p and (13. ip . Therefore, we finally obtain 



m 2 3m (A - Sm^) . 5 {?,A^m - 2SArn? + 63m^) 

^ 128y4iV2 



105 (m (A^ - 21A2m2 + 99Am^ - 143m6)) „ 

^ ^ ^A^ + --- (3.10) 



2048^15/2 

Note that ai vanishes in the limit of A — )■ oo so that a ~ \/~A asymptotically. This is 
natural because A = m^ + C2 — ^ is identified with the VEV of the adjoint scalar [2T] 
and the limit A oo corresponds to the UV limit in gauge theory. 
The free energy (12. 3p can be obtained using the relation 

^^^matrix ^ v^A / A, + ^ ) = + 2^7 v^A(^ A,) . (3.11) 



II Similar structure has been noticed in 1301 in the semi-classical approach of the system for each 
iVi^'s. 
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Noting the asymptotic behavior of the resolvent W{z) ~ ^/PgN/z + \/^g{J2i '^i) /-^^ + 
0{l/z^) one finds the terms of order 0{z~'^) in the loop equation fl2.12p 

2g^P{Y,h) = C2 + {fil-fil) + 0{g'). (3.12) 
Since we put fii = fi^ = m, fl2.12p reduces to 

^M^""*'" = 2c2 + C(/) = 2{A-m^)+A^ + 0{g^). (3.13) 

To integrate the right-hand side with respect to A, one needs to find the explicit 
form of A. Using the filling fraction a = Oq + ^cli given in (13. 7p and (I3.10p one finds 

^ = H -A^ H ^ A^ 

' 2a2 32a6 64aio 

Sa^ - 252a^m2 + 1638aW - 2860a2m6 + 1469mS , „ 

A^ + ■ ■ 



14 



8192a 

e+ / fn ^r, Sa^m — 5m^ . . 5a^m — 280^771^ + 27m^ . « 
H — A H A A^ 



6 32aio 



63a^m - 819a^m^ + 2U5a'^m^ - 1469m^ . „ 

+ Mi^' ^ + 

+0{/). (3.14) 

Putting this result into fl3.13p . we obtain the free energy of the form i^™a,trix _ ^matrix _|_ 
^^matrix _^ 0{g'^) wherc Fq™''*"'' is the planar contribution 

Fo'"'^*""(a,m) = 2(a^-m^)logA+ A^ + A^ 

5a^m^ — 14a^m^ + 9m^ . « 

H A^ 

192aio 

5a« - 252a^m^ + 1638a^m^ - 2S60a'^m'^ + U69m^ . „ 

H A*^ H 

32768ai4 

(3.15) 

and F^'^^'^^^ is the half-genus correction 

™ntri^/ X "^.9 3a^m — 5m^ . 4 Sa'^m — 28a^m^ + 27m^ . 

^matrixj^^^ m) = -A^ H A^ A^ 



a 



16a6 96aio 

63a6m - 819a^m3 + 2145a2m^ - 1469m^ . „ , , 

^ 77^?^. ^ +■■■■ (3-16) 



4096a 
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Here F^'^^^^^ is the same as evaluated in [21] and agrees with the Seiberg-Witten 
prepotential of the corresponding gauge theory, while F^^^^^^ perfectly agrees with the 
half-genus correction to the Nekrasov partition function with parameters /ii = /is = m. 
(The corresponding Nekrasov partition function is given in (IC.lSp for comparison.) 

4 Np = 3 case 

For Np = 3 case, the corresponding potential of the matrix model is given by (12. 8p 
with the definition (l2Aill 

Vo{z) = jj^slogz + mi\og{z -1) - Vi{z)=\ogz. (4.1) 



f{z) has the following expression 



m = ^ + ^ + % (4.2) 

z z — 1 z^ 



and Cl , C2 and C3 are given by 



^ / A \ 

The symmetry of the partition function V^'(A/)) = leads to the condition 

Cl + C2 = 0. (4.4) 

Another constraint on C2 and C3 follows from the asymptotic behavior of the loop 
equation (I2.12p where the resolvent has the asymptotic behavior Wi^z) = y/PgN/ z + 
0{1/ z'^). Noting that the derivative of the potential is V\z) = (7x3 + rrii + e+/2)/z + 
0(1/ z'^) and f{z) = (c2 + Ca)/^^ + 0{l/z^) with ci + C2 = one has the left-hand side 
of the loop equation at the order of O^l/z"^) 

C2 + C3 = 4:^/PgN [^/l3gN + /is + mi) = - (/is + niif , (4.5) 

where in the second equality we used equation (12.71) . The two constraints (14. 4 p and 
(14.50 allow one independent parameter, for which we choose C3. 
The planar spectral curve has the form 
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where P4,{z) is a polynomial of degree four which is given by 

P^iz) = 4m^/ - 4(5 - miA - + m^)z^ 

+ (45 - 4miA + - 4/i3A)z2 + 2A(2/i3 - A)^ + A^, (4.7) 

with B = C3 — ^3 A + yUg. Hereafter, for simplicity, we set mi = rrioo = and /X3 = m 
as in [21j. Then 54(2;) becomes a third order polynomial of z as 

P4(^) = -45;z^ + (45 + A^ - 4mA)22 + 2A(2m - A)^ + A^ (4.8) 

The Coulomb branch parameter identified as (I2.24p has the planar contribution 



ao 



/ m . B + 3m^ 



6452 



A2 + 



m 



25653 



(Sm^ + 5)A^ 



1 



(35^ + 30m^5 + 175m^)A^ 



1638454 
m 



655365 



-(952 + 70m25 + 441m^)A^ + 0{A^) 



(4.9) 



and the half-genus correction ai (12.271) has the form 

A / 1 /I 
Oi = (b — , 1 ] dz. 

27r« J A Jl\(z) \z 



As in Np = 2 case, one has 



by noting 



ai 



dan 



dm 



+ A 



dan 



B 



dB 



dm 
daQ 



dB 



B: fixed 



m: fixed 



A /■ dz 
27rz J A z.fPlz) 



1 



dz 



(4.10) 



(4.11) 



(4.12) 
(4.13) 



where we treated m and 5 in (14. 8 p as independent variables. Hence, ai is given as 
follows: 



ai 



A ^ mA2 (5 + 3m2) A^ ^ (95m + 25m3) A^ 



A^/B 3253/2 25655/2 



(952 + 905m2 + 245m4) 
6553659/2 



A^ + 



m 



409657/2 
(7552 + 4905m2 + 1323m4) 



524288511/2 



A^ + - 



(4.14) 
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Note that Oi vanishes in the hmit 5 — )■ oo, which guarantees the asymptotic behavior 
a ~ oq. 

Now, we evaluate the free energy of the matrix model 

A-^F"^"'"" = -2v/^^V/^\ = C3 = B + mA-m\ (4.15) 



dA 

In order to integrate this, we solve the inverse function of a{B,m) as 

„ , 9 niA + . n a'^ — 6m^a^ + 5m^ . ^ 5a^m^ — 14a^m^ + 9m^ . ^ 

B = { \ —A"' H A^ H A^ + 

' 32a2 8192a6 262144aio 

e+ / A ni , r, ^a^m — 5m^ . . 5a^m — 280^771^ + 27m^ . ^ 

+— H — A"* A^ + 

2^2 16a2 2048a6 131072aio 

+0{g'). (4.16) 
By integrating f l4.15p one obtain the free energy of the form 

— Gm^a^ + 5m^ . . 5a^m2 — Ma^m^ + 9m^ . ^ / _x 

+ ^T^;;7:7r^^ + .r..ar. in ^ +■■■ (4-17) 



32768 a6 1572864 a^o 

^^„tH^ A Tn ^r, Sa'^m — 5m^ . . 5a'^m — 28a^m^ + 27m^ . f. 

pmatnx ^ ^^2 _^ ^4 ^6 _^ 

^ ^ 32a2 8192 a6 786432 a^o 



(4.18) 

Here i^™^*"'^ was evaluated in [2T] and is equal to the Seiberg-Witten prepotential of 
Np = 3 gauge theory. On the other hand, the half-genus correction i^™a,trix newly 
evaluated here and perfectly coincides with the half-genus correction to the Nekrasov 
partition function of SU{2), Np = 3 gauge theory for our parameters mi = moo = 
0, yUs = m (See flC.131) in appendix [Cl) . 

5 Np = 4 case 

For Np = 4 gauge theory, the matrix model is given by (12.51) with 

Vo(-2) = mologz + mi log(2 - 1) + m2 log(z - g), 1^1(2;) = log 2;, (5.1) 

so that V{z) = Vo{z) + {e+/2)Vi{z). The function f{z) defined in (I2.13P is now written 
as f{z) = XlLo where q's are given by 



Co = _ 

12 



From the symmetry of the system '^j{V'{\i)) = 0, it follows that 

2 

^Q = 0. (5.3) 

2=0 

From the residues at infinity, we find an another constraint 

Ci + qc2 = 4:^/^gN (^/^gN + =mlo- ' (^-4) 



i=0 / \i=0 



where in the second equality (12. 7p is used. Thus, essentially there is a single free 
parameter, for which we take Cq. 

The planar spectral curve is exactly what was obtained in 



where P4,{z) is a degree four polynomial of ^j**] Hereafter, setting mo = rrioo = 
and mi = m2 = m for simplicity (which provides the equal mass for hypermultiplets 
/ij = m for i = 1, ■ ■ ■ , 4), one has P^iz) of the third order polynomial 

P^z) = Cz^ + {{1-qfm^ -C{l + q)} z^ + Cqz, (5.6) 

where we defined C = qcQ. 

The Coulomb branch parameter is given by (I2.25P . The planar contribution was 
already calculated in [21], which is written in our notation as 

, /— / h2{q)m^ h^{q)m^ h^{q) [ m"^ 

Co = iyC[hQ[q)-hi[q)— t^t^ ^7^7 + C ' 



C 3 C2 5 C3 7 C4 ' V C5 



(5.7) 



where hi{q) is defined in terms of the expansion coefficients of a hypergeometric func- 
tion ^ 

"^('l - ^i'^-ds^'^ 131572'' ("°) 

I. / N 5 125 1125 2 125 , 125 4 ^^,5, 

/i3(g) = \ qH q^ + Oiqn. (5.11) 

16 64 ^ 1024^ 4096^ 262144^ ' ^ ' 



*For the explicit expression of Pi{z) for general mass parameters, see equation (3.39) in [21) . 
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The half-genus correction oi is given by fl2.27p 
m(l + q) f dz m(l — g)^ 



This integral is expressed in terms of ag: 



dz 



{z - l){z - q) 



. (5.12) 



dan 



dC 
daQ 



dz 



dm 

Hence, we finally obtain 



m: fixed 47ri Xl 

m(l — f dz 

C: fixed 



27ri 



^/p^{z-l){z-qy 



(5.13) 
(5.14) 



ai 



2m(l + q) 



dan 



dC 



dan 



dm 



c 



^ 7 
Tfi T71 771 771 

+ 93{q)-^ + qM-^ + ariq)-^ + 



(5.15) 



where gi{q) are functions of q defined in terms of hi{q). The first few components are 
given by 



9i{q) = {-ho{q) + 2hi{q) - ho{q)q), 
Ah2{q) - 3{1 + q)h{q) 

g-Aq) -- 
9b{q) -- 
97{q) -- 





3 


Qh{q) 


-5(l + g)/i2(g) 




5 


8h^{q) 


-7{1 + q)h{q) 
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(5.16) 
(5.17) 
(5.18) 
(5.19) 



To find the free energy F'^^^'^^^ of the matrix model one can use the equation 
d 



^matrix _ 4g^^^ Y^^^ ) = 4:m2W{z) 



\z=q 



From the loop equation (12.121) . we find that W{z) in the vicinity of ^ = g 

C2 



W{z) 



4ms 



+ Oiz-q), 



(5.20) 



(5.21) 



which implies that 



d_ 
dq 



nmatrix 



C2 



:i-q) 



Am' 



c 



(5.22) 
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Here C is obtained from a{C) as 
C{a) = -0? 



1 2hi{q)m^ 2ho{q)h2iq) - 3hi{q)^ m 
ho{qy hoi^q) 3 a 



4 



lQho{q)hi{qf - 10ho{q)^h{q)h2{q) + 2ho{qfh3{q) 



m / 2(l + g)/io(g)-4/ii(g) 8/io(g)/i2(g) - 12/ii(g)^ 



^o(g) 



+ 



12/io(g){5/ii(g)3 - 5/io(g)/ii(g)/i2(g) + /^o(g)^/^3(9)} 

5 

This is used to integrate f l5.22p to get the free energy _pmatrix 



+ 



(5.23) 



Tjimatrix 
^0 



(a^ — m^) logg + 



2a2 



(13a8 + lOOa^m^ + 22a^m^ - 12a2m6 + 5m^) ^ 
^ 64^6 ^ 



23ai2 ^ 204ai°m2 + ^la^m'^ - A%a^m^ + 45aW - 2%a^m^^ + gm^^ 



192aio 



32768a 



14 



(2701a^^ + 2UA<da^W + 71QAa^V - 9000a^°m^ 



^1 



matrix 



+ 12190a*m* - 13384a^m^° + 10908a^m^2 - 5720a2m^^ + 1469m^^) 
+0{q'), (5.24) 
2m(a2 + m^) 9a^m + lla^'m^ — ga^m^ + 5m'' r, 

—o -Q TT-. (1 



SSa^^m + 51a^m3 - 72a''m5 + 90a^m^ - 70a^m^ + 27m^^ , 

q 



48a 



10 



2048ai4 



(I257a^^m + 1791a^^m^ - 3375a^°m^ + 6095a®m^ - 8365a6m^ 
+8181aW^ - 5005a^m^^ + 1469m^^) q^ + C(g^). (5.25) 



Here i^matrix ^g^g calculated in [21j and shown to be equivalent to the Seiberg-Witten 
prepotential of Np = 4 superconformal gauge theory. The newly found half-genus cor- 
rection i^^^ti'i^ coincides with the half-genus correction to the corresponding Nekrasov 
partition function fIC.lip for /ij = m. 
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6 Summary and Discussions 



In this paper, we investigate the /3-ensemble of the matrix model known as Penner type 
matrix modeL By solving the loop equation generalized to /3 7^ 1 we explicitly evaluate 
^matrix rpj^^ result perfectly agrees with the Nekrasov partition function with general 
i7-background parameters €1,62, corresponding to the M = 2,SU{2) gauge theories 
with Np = 2,3 and 4 flavors. 

It is noted that the relation of the filling fraction with the Coulomb branch param- 
eter is to be modified as in (12.241) at the order of half-genus expansion, which reduces 
to the original one proposed in [13j in the planar limit. The modification is understood 
as a condition of the UV limit of the gauge theory so that the the correction to the 
Coulomb branch parameter need to vanish at each order of (half-) genus expansion as 
the corresponding scale parameter becomes infinite. 

The loop equation provides a systematic way of finding the (half-) genus expansion 
of the theory. The merit of the expansion is that even at this first non-trivial order, 
the higher instanton contribution can be obtained by the simple algebraic algorithm. 
To go to high genus expansion, one needs to solve the lower genus result of multi-point 
resolvent and the result will be reported in the near future. In addition, the extension 
of the analysis to the multi-matrix model is highly desired so that the Selberg integral 
approach of the SU{N) gauge theory [12] is to be compared. 
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A AGT conjecture and Matrix model 

We briefly review the AGT conjecture [1] and its relation with the Penner type matrix 
models which is viewed as dual to d = 4, J\f = 2 gauge theories [I3j for the setup of 
our convention in the manuscript. 

A.l AGT relation 

The authors of [1] pointed out that Nekrasov partition functions of d = 4,A^ = 
2, SU{2) gauge theories is identified with conformal block of of the four-point cor- 
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relation function of the Liouville theory, 



Liouville 



c(/x/,ei) 



2 

q 2 ^0+^ 



a da 



Vkrasov(«'/^/'ei) (A.l) 



where = e^"'^^^-' is a vertex operator of dimension Aq, = a{Q — a). Q = b + 1/b 

is the Liouville background charge and b is the parameter in the Liouville potential 
e^^*^. -^Nekrasov 5 /^^' ^i) right-hand side is the Nekrasov partition function of 

A/" = 2, Np = 4, SU{2) gauge theory and depends on the Coulomb branch parameter 
a, masses of four hyper multiplets fij, and the ^-background parameters ei and €2- 
The prefactor c(/i/, e^) depends only on /i/ and e^. 

g in flA.ip is the exponential of the UV gauge coupling and the relation (lA.ip is 
complete if the Liouville parameters are given in terms of the gauge field parameters. 
If one introduces a mass scale h (which was set to be 1 in [T]) so that hrha = ma, one 
has 

Hi=mi+moo, iJ,2 = mi-m^, /is = m2 + mo, /i4 = "^2 - mo . (A. 2) 

f2-background parameters are identified with the Liouville parameter b as 

= hb, t2 = h/b, (A.3) 
so that the background charge is written as Q = (ei + £2)//?- = e+/^- 

A. 2 Dijkgraaf-Vafa's proposal of matrix model 

Inspired by the AGT relation, a matrix model description was proposed in [13]. Note 
that the Liouville correlation can be evaluated perturbatively 

flj dXjdXj e'"^^''^^ V^^^<,{oo)V^,il)V^MVf^^^^iO)j (A.4) 

where the expectation values is evaluated in terms of the free field description {(j){z)(f){w)) 
-l/2\og{z-wY OT (e2"i</'We2"2 0(«')^ = \z-w\-^°'^°'\ The number of integrals comes 
from the perturbation of the Liouville potential and is viewed as screening integrals. 
The non-vanishing contribution is obtained if the number of screening charge satisfies 
the neutrality condition 

2 

^m + m^ + bN = 0. (A.5) 

1=0 
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The correlation ( ]A.4p is given as (up to the g-independent prefactor) 



2(mo+<;^/2)m2 : 

q (1 _ q y 



where 



N 



.1=1 



KJ 



2b ^ 



V{z) 



mQ + — ]\ogz + mi \og{z - 1) + m2 log(z - q). 



(A.6) 



(A.7) 



Note that /4 is not well defined unless the integration ranges and parameters b and 
h are to be appropriately arranged to make I4 convergent. To fix this problem, one 
may consider the integrals with 6^ = — and hb = 2g^/P (or b = h = —2ig) so 

that integration is well defined even when any two of the integration variables coincide. 
This defines the partition function of matrix model 

N 



^matrix 



.1=1 



KJ 



-A., 



exp 




(A.8) 



where A/ is the eigenvalue of the hermitian matrix and the size of the matrix is 
given from the neutrality condition (1A.5I1 or (12. 7p . However, the matrix model is not 
the usual one unless /3 = 1 and is called /3-deformed matrix model, or Penner type 
matrix model. 

As seen in section |2l ^matrix depends on a single parameter which is not in the 
potential V{z) and the parameter can be chosen as the filling fraction. The origin 
of the ambiguity comes from I4 in (]A.6|) where one needs to arrange the integration 
range appropriately, whose origin also traces back to the perturbation expansion of 
the Liouville correlation f lA.4p . This ambiguity does not appear if one evaluate the 
Liouville correlation using the conformal block as put in ( lA.ll) . To resolve this discrep- 
ancy, one may view the perturbation as the one with the fixed filling fraction so that 
the integration ranges are chosen so that Ni number of integration ranges from to g 
and N — Ni number of integration ranges from 1 to 00 resulting the filling fraction is 
related with the Ni/N. 

According to AGT, VEV of SU{2) gauge group is identified with the Liouville 
momentum a of the conformal block in ( IA.1|) . On the other hand, the Seiberg-Witten 
curve is "quantized" and the VEV is identified with the filling fraction of the S-W 
curve. In this sense, it is very natural [13] that the Coulomb branch parameter of 
the gauge theory is identified with the filling fraction of the matrix model. With this 
interpretation in mind, one may identify 2'matrix with ^Nekrasov as follows: 



(mQ-m2) -2m2e^ 

q ^ (1 -g)" 



7SU(2) 
^Nekrasov" 



(A.9) 
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On the other hand, it was pointed out in [T] that the Nekrasov partition functions 
of SU{2) and U (2) gauge theories are related by 

r7SU(2) I . (Ml+M2)(M3+M4) TJ ( K ^ r\\ 

VkrasovK/^/'ei) = ^ Vkrasov(«'/^/'ei), (A.IO) 

Note that the f/(2) gauge theory has two independent Coulomb branch parameters 
a = (ai, 02), but we set a\ = —a2 = a in the right-hand side of ( lA.lOp . Note also that 
our /i/ are masses of anti-fundamental hyper multiplets, as states in appendix O By 
combining (1A.9P and (]A.10|) . we obtain 



^matrix = Q ^Nckrasov' (A-HJ 



up to a prefactor which is independent of q and a. 



B /5-deformed version of loop equations 



We present the derivation of the loop equation for /3-deformed Penner type models, 
following [28] (See also [261 129]). We start from the partition function of the form 

N 



l[dXj 



1=1 



(B.l) 



and consider the change of the integration variable as A/ — > A/ + This changes 

the expression of the integrand as well as the measure. Collecting terms proportional 
to e, we obtain 

/ ^ K N 







E 

/,j=i 



1=1 



1-/3 



(A7 



TV 



1 



9 ''i^,\z-h 



/V\z)-V\h) 

9 h 



z — \i 



(B.2) 



Here, the first and second terms come from the variations of the measure and A^, 
while the third and fourth terms are from the variation of the potential V^(A/). By 
defining multi-point (connected) resolvent as 

2-s 



W{z,, 



^Zs)=P 



^1 — A 



1 



2, — A/ 



(B.3) 



h ' Is 

the equation ( IB. 21) gives the loop equation for the /3-deformed matrix model 



= g^W{z,z) + W{z)' + g(^,/p-^^W'{z) + V'{z)W{z)-^{BA) 



N 



1=1 



V'jz) - V'jXi) 
z — \j 



(B.5) 
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Note that the first term in (16.40 is 0{g'^)^ while the third term is 0{g) which con- 
tributes to the half-genus correction. Omitting higher order term, one has the loop 
equation up to half-genus correction with 6+ = 2g{^/P — 1/v^) 



W{zf + ^-W'{z) + V'{z)W{z) - ^ = 0. (B.6) 



C Nekrasov Partition Function 

We present the Nekrasov partition function [2], [3] of f/ (2) gauge theories for the com- 
parison with the matrix result. The Nekrasov partition function is factorized into the 
following three contributions: 



^ class 



^1 — loop -^inst ) 

(C.l) 



where Zdass and Zi_ioop are the classical and 1-loop contributions, respectively. In the 
following, the instanton part is elaborated for Np = 4, 3, 2 cases. 



C.l Nf = 4: theory 

(C.2) 

where the sum runs over pairs of Young diagrams Y = {Yi, Y2) and a = (ai, 02) denotes 
the Coulomb branch parameter while fii are mass parameters of four fundamentals. 
The vector multiplet contribution Z^cc{S,Y) is given by 

2 

Z,ecia,Y) = II II E{a^ - aj,Yi,Yj, s) l[{e+ - E{a, - a^,Yj,Yi,t)), {C.3) 

i,j=l s(^Yi teYj 

where s and t run over all the boxes in Y^ and Yj, respectively. The constituent 
E{a, Yi,Yj, s) is defined by 

E{a,Yi,Yj,s) = a-eiLY^{s) + e2{AYAs) + l). (C.4) 

The arm-length Ay{s) and leg-length Ly{s) are defined by 

AY{s) = \k-l, LY{s) = \'i-k, (C.5) 

where {k, I) is the coordinate of the box s in Y. The two integers and A'^ denote 
the heights of k-th column of Y and Y'^ , respectively. 
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The anti-fundamental hyper multiplet contribution is replaced by 

2 

^afund(a,>^,yu) = JJ + (C.6) 

where 0(a, s) = a + ei{k — 1) + €2(1 — 1) when s is on the position (fc, I) in Y. Note 
that in |2T] the last two factors in (IC.2P were given as 

where Zfy^^^{a,Y , ^) is a contribution from a fundamental hypermultiplet with mass 
/i, which is equivalent to our expression if e+ = 

^afund(a, 'i^, /i3)^afund(a, /i4) (C.8) 

— * — * 

because Zafund(o, yu) = ^fund(a, ^, e+ - Z^)- 

From now on we concentrate on the case of a = (a, —a) of interest. Defining the 
free energy of instantons by 

F^-^* = (-eie2)logZinst. (C.9) 
we can expand _F™*^* in powers of g (with the relation (12. 2p ) 

In addition, if we put /i/ = m as treated in the text, we have 

^" = « + ei^^ « 



23ai2 204ai°m2 + 51a*m^ - 48a6m6 + 45a^m* - 28a2mi° + Om^^ ^ 

192aio ^ 

1 



(2701a^^ + 26440a^^m^ + 7164a^^m^ - 9000a^°m^ 



32768ai4 

+12190a*m* - 13384a^m^° + 10908a^m^2 - 57200^771^^ + 1469m^^) 



jj^g^. 2m (a^ + m^) 9a^m + lla'^m^ — 9a^m^ + 5m''' 2 
= a- ^ 8^^ ^ 



38a^°m + 51a*m3 - l^a^nt' + 90a^m^ - l^c?m^ + 27m" 



10 ^ 



48a 

^ (I257a^^m + 1791a^2m^ - 3375a^°m^ + 6095a*m^ - 8365a^m^ 



2048ai4 

+8181aW^ - 5005aW^ + 1469m^^) + C(g^). (C.ll) 
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C.2 Nf = S theory 

The instanton part is given by 

A3 'Zvec(a, F)^afund(a, /^l ) ^afund («, Y, /i2)^afund(a, F, /i3),(C.12) 

y 

where /ii,/i2 and /is are the masses of three hyper multiplets and A3 is a dynamical 
scale. The Coulomb branch parameter a generally has two independent components. 
However, as we have seen in appendix El we only need to consider the case of a = 
(a, —a) in this paper. The half-genus expansion of the free energy has the form 

pmst ^ _^ _^ A 2 _^ ^ ^4 ^ ^ ^6 _^ 



2 64a2 32768a6 1572864aio 

pinst = -A^ H ^ — ^A^ ^ -A^ + 

^ 2 32a2 8192 a® 786432 a^o 



(C.13) 



where we set /ii = /i2 = and /is = m. Note that the first term in Fq and Fi is 
independent of the Coulomb branch parameter a. 

C.3 Nf = 2 theory 

The instanton partition function is given (with a = (a, —a)) 

^inst = ^ Aa'^'Z^ecla, l^)^afund(a, ^, /il)^afund(a, Y, /is). (C.14) 

y 

The half-genus expansion of the free energy has the form 



-|- . r, a'^ — Ga^m^ -|- Sm"^ . ^ ba^m? — lAa'^m'^ + 9m® 

A H A H 

2a2 64a6 192 a^o 



5a^ - 252a®m2 + 1638a^m^ - 2860a2m'' + 1469mS . „ 

H A + ■ ■ 

32768 



^inst _ _-/^2^::is:z — _ ^^^^ — ' "'"^ ^ A 



m 2 m (3a2 — Sm^) . m (5a^ — 28a2m2 + 27m 



4^ 



16a6 96aio 



m (63a® - 819a^m2 + 2145aW - 1469m®; g 
+ 4096^^^ ^ ^^-^^^ 



where we set /ii = /is = m. 
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